Generalized Bose-Einstein phase transition in large-m component spin glasses 
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It is proposed to understand finite dimensional spin glasses using a 1/m expansion, where m is 
the number of spin components. It is shown that this approach predicts a replica symmetric state 
in finite dimensions. The point about which the expansion is made, the infinite-m limit, has been 
studied in the mean-field limit in detail and has a very unusual phase transition, rather similar to 
a Bose-Einstein phase transition but with N 2 ^ 5 macroscopically occupied low-lying states. 
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After almost three decades of research, the nature of 
the low temperature phase of finite dimensional spin 
glasses is not understood. The usual approach is to start 
from the exactly soluble Sherrington-Kirkpatrick (SK) 
mean field model and expand about it towards finite 
dimensions. This method resulted in the monumental 
replica field theory, summarized in 0] by three of its main 
contributors. The results of this theory are that replica 
symmetry is broken in finite dimensions, just as in the 
SK limit, down to 6 dimensions, below which all known 
calculational tools break down and not much is known 
analytically. Yet there exists a mathematical proof that 
this picture cannot hold in any finite dimension Q . This 
proof is unfortunately non-constructive, so it does not 
give any insight into the nature of the spin glass phase 
and it cannot decide between alternative scenarios such 
as the droplet or the chaotic pairs picture. Neither does 
it show where the replica field theory goes wrong. We 
therefore propose an alternative method to investigate 
the finite dimensional spin glass phase in order to by- 
pass conventional replica field theory and to obtain a 
theory which does not contradict the exact mathemati- 
cal results. To this end we shall expand about the infi- 
nite component limit of the m-component spin glass in a 
power series in 1/m. As we will see, an expansion of this 
type preserves the replica symmetry found in the m = oo 
model and will therefore give a picture of the spin glass 
phase which is quite different to that of the usual replica 
field theory. As part of our programme it is necessary 
to analyse the m = oo-component SK spin glass phase. 
The low temperature phase of this model is a generalisa- 
tion of a Bose-Einstein condensation in the sense that the 
spins condense into a no-dimensional subspace of the Tri- 
dimensional space they can occupy, where no is a number 
of order TV 2 / 5 @. 

We study the spin glass model defined by the Hamil- 
tonian 



H 



- JijSi.Sj hi 



(1) 



where the N spins Sj are vectors with m components and 
we use the normalization s 2 = m. The vectors h; are m- 



component gaussian random fields with correlator 



hfh) = h 2 dijd ab 



(2) 



and field strength h. 

In this paper our numerical work is on the mean-field 
limit where the off-diagonal are independent gaussian 
random variables with variance l/N, and J a — 0. Some 
of our results, however, also apply for corresponding 
to a finite dimensional system. In all cases the partition 
function at a temperature T — 1/(3 can be written as 



Z = 



0H 



(3) 



-in 



x exp I -^y^Hj(m- sf) + - J^Sj-Sj + h^.s. 
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where we have introduced an integral representation for 
the S functions. The integrals over the spin components 
sf can be done, and making use of Eq. J2J this results in 



Z = 



x exp 



n 



dHi(3 
An 



where the matrix \ is defined by 

Xij = {A^ 1 )^ with 



, (5) 



(6) 
(7) 



For large m, the integral in Eq. (JjJJ can be evaluated by a 
steepest descent calculation in the Hi, giving rise to the 
conditions 



(3 = Xu+Ph 2 (x 2 )i 



1. 



(8) 
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which determine the Hi. Eq. JSJ can be rewritten using 
the eigenvector decomposition of A, the inverse of Xi as 



E 



(a 



n\2 



A, 



/to 2 



(9) 



where X n are the eigenvalues and a™ are the orthonormal 
eigenvectors of A. 

At zero temperature, Eqs. JHJ and JjJJ are no longer 
well-defined. In the ground state all spins are aligned 
parallel to their local field, i.e. 



HiSi 



(10) 



For the large-m limit there exists a unique stable solution 
of these equations for the Hi, as opposed to the case of 
the finite m spin glass which has an exponentially large 
number of stable solutions Q . This fact allows us to solve 
Eqs. H10|l numerically without difficulty We note that it 
is not immediately obvious that the Hi in Eq. I|10fl are in 
any way related to the Hi buried in the matrix x i n the 
finite temperature problem, Eq. ©. It has been shown in 
, however, that the Hi in Eq. IjlOl) are equal to the limit 
of the Hi in Eq. JHJ as (3 — > oo . This observation allows us 
to cover the complete temperature range including T = 
within one framework. 

Given the Hi determined by Eq. © and disregarding 
any irrelevant prefactors, the partition function is then 



^ ^2(H t + h\ lt ) + i lndet(x//^ 
x exp (-^\ndct(B/[3 2 ) - — mm) , (11) 



exp 



where we have included the first order fluctuation cor- 
rections around the steepest descent values for Hi which 
involve the matrix B defined by 



Bi 



1 



(12) 



This expansion is valid for any spatial dimension d in- 
cluding oo since we have so far made no assumptions 
about the matrix ,7^. 

It has been shown by de Almeida et al. that the 
to = oo SK spin glass is replica symmetric and in the 
thermodynamic limit has the same free energy as the 
spherical spin glass model Q (but despite having the 
same free energy, the physics of the low temperature 
phase is in fact very different from the spherical model). 
All the Hi are equal in the thermodynamic limit, and the 
distribution of eigenvalues of the matrix A follows the 
Wigner semicircle law. The phase transition (at h = 0) 
then follows from the fact that Eq. JHJ only has a solution 
for < (3 < P c = 1, which is 



Hi = {3 + l/[3. 



(13) 



At the critical temperature T c = 1/(3 C > the smallest eigen- 
value becomes zero. 

If the field is non-zero, however, there is no phase tran- 
sition since the term involving 1/A 2 allows for a solution 
at any temperature. This can be demonstrated by eval- 
uating Eq. © under the assumption that again all Hi 
are equal to H in the thermodynamic limit and thus the 
Wigner semicircle law holds, which then yields after some 
algebra 

p 2 h 4 = (p 2 + 1 + [3 2 h 2 - f3H)(H 2 -4). (14) 

This equation has a positive physically relevant solution 
H for any (3, as long as h ^ 0. 

The scenario in finite dimensions has been described 
in 0- This work was an extension of the Bose glass 
theory of Hertz, Fleishman, and Anderson 8]. For tem- 
peratures T > T c the eigenfunctions corresponding to the 
smallest eigenvalues of the matrix A are localized but be- 
come extended at criticality. The Hi vary from site to 
site. However, Eq. ijHJ) still applies. The second term 
on the right hand side of Eq. JSJl involves the replicon 
susceptibility x 2 ■ Assuming there is a phase transition 
in non-zero field.this is the quantity which diverges on 
the critical line |9j. This, however, makes it impossible 
for Eq. ijHJ) to be satisfied, thereby ruling out the exis- 
tence of a phase transition in a field, i.e. an Almeida- 
Thouless line ^(J; by contradiction. As this line marks 
the onset of replica symmetry breaking, we deduce from 
its absence that the large m limit and the straightforward 
expansion in 1/m about it, produces a theory which is 
replica symmetric. The absence of such a line is a pre- 
diction of the droplet theory of spin glasses. However, 
usually with 1 /to expansion methods one does not adopt 
the direct expansion approach in powers of 1/m but in- 
stead uses (ad-hoc) self-consistent approximations, such 
as could be generated by making the sum of the terms in 
to and of O(l) in the exponentials in Eq. Qllf) stationary 
with respect to the Hi, and it is possible that with such 
an approach replica symmetry breaking might emerge. 

We turn now to a more detailed study of the phase 
transition mechanism at zero field in the SK limit. First 
we examine the ground state properties in zero field. 
Eq. I|10fl can be written as 



3 



{H i 5 ij -J ij )s?=Q ) 



(15) 



where a = 1, 



, to labels the spin components, show- 



ing that the matrix Ajj = HiSij — Jy has at least one 
eigenvalue which is exactly 0. On the other hand, it was 
shown in Q that there is an upper limit of the number 
of null eigenvalues no, which is no < y/2N. 

The number of null eigenvalues will turn out to be the 
key quantity of the low temperature phase of the large-m 
spin glass. Here, we can already see some implications. 
The fact that there are no null eigenvalues reduces the 



3 



effective number of spin components from m to no. To 
see this, consider the matrix formed by the entries of the 
spins, sf. Since the rows of this matrix (regarding a as 
the row index) correspond to null eigenvectors of A and 
there are only no linearly independent ones, the row rank 
of this matrix is at most no. But since the row rank and 
the column rank of any matrix are equal, the column rank 
is also at most no. Therefore the N columns, being the 
spins, can only span a no-dimensional subspace of the m 
dimensional space they live in. Making use of the global 
rotation invariance of the spins, it is therefore sufficient 
to set to = no, or at least m = \/2N > uq when no 
is yet unknown. Obviously, this accelerates numerical 
simulations, allowing one to go to relatively large system 
sizes, as we will see in the following. 

We have solved Eqs. i|10|) numerically for the spins s, by 
straightforward iteration, setting Hi — | Y^,j Jij s j\/V™ 
at each step, until the average angular deviation of the 
spins from their local field directions was smaller than 
some prescribed accuracy (we used 10~ 8 ). From the 
ground state configuration found in this way we deter- 
mined the matrix A and analysed its eigenvalues. We 
found that the smallest non-null eigenvalues are always 
at least 5 orders of magnitude larger than the null eigen- 
values (which we find numerically to be of the order 10 -8 
or smaller). The average number of null eigenvalues as a 
function of the number of spins is plotted in Fig. In 
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FIG. 1: The average number of null eigenvalues (no) as a 
function of system size TV. The error bars on the data points 
are smaller than the point size. The dashed line is ~ N 2/5 for 
comparison. 

the range of system sizes accessible to us, it behaves as 
no ~ TV 2 / 5 . This agrees with the prediction from Q . The 
calculation leading to this result was, however, based on 
an approximation of the density of states of the matrix 
A which, at finite values of N, had a gap in the den- 
sity of states and a square root singularity at the edge of 
the gap. Numerically we found a very different behavior, 
see Fig. [21 Instead of a gap there is an enhancement in 
the density of states at low eigenvalues (rather similar 




FIG. 2: Density of eigenvalues of the matrix A for a system 
with N = 450 (but typical for all system sizes) , excluding the 
null eigenvalues, averaged over 150 samples. The density p(A) 
goes to a constant at small A. (The constant tends to zero 
for large system sizes, and for an infinite system the Wigner 
semicircle is restored.) 



to the behavior produced in the three-dimensional XY 
spin-glass model by finite size effects 11]). 

A new argument is therefore needed to explain the 
observed value of the exponent 2/5. From 0, |j| it 
can be deduced that the ground state energy per spin 
of the m-component spin glass, divided by to, goes as 
— 1 + l/4m + C(1/to 2 ). Since no is equal to the effec- 
tive number of spin components, the system would at 
first sight be able to attain its lowest energy state by 
choosing no as large as possible, i.e. equal to its upper 
bound. However, this calculation of the energy was done 
taking the thermodynamic limit first, such that m (or 
n ) is always much less than any power of N. When 
no is comparable to some power of N, there are addi- 
tional energy costs, whose magnitude can be estimated 
by the following argument. Starting from all Hi being 
equal and then tuning them in such a way that there 
are no null eigenvalues will result in a downward shift in 
the eigenvalue spectrum of A which is of order (no/-/V) 2 / 3 
(the lowest no eigenvalues in a Wigner semicircle reach 
this far from the band edge, and the shift is expected 
to be of the same order). The ground state is reached 
when the energy arising from these two competing terms 
l/4no + const. (no/iV) 2//3 is minimized, which is the case 
when n ~ TV 2 / 5 . 

It is striking that in the ground state the individual 
spins condense into a 7i -dimensional subspace of their 
original m-dimensional space. This behaviour is a gen- 
eralization of the conventional Bose-Einstein condensa- 
tion where the constituents condense into a single (one- 
dimensional) state, which is also what happens in the 
spherical spin glass model Q. 

To show that the behavior we observed at zero tem- 
perature is in fact the generic behaviour in the low tem- 
perature phase, we solved Eqs. J8J numerically in zero 



field. Since there is no phase transition for a finite num- 
ber of spins, we were able to solve Eqs. JSJ numerically 
over a very large temperature range, both in the high and 
low temperature regions. We used a standard Newton- 
Raphson iteration scheme for this purpose. This method 
converges very quickly for high temperatures j3 < 1 and 
fails to converge for low temperatures unless the starting 
configuration for the iteration is already sufficiently close 
to the solution. In order to ensure this, we employed the 
exact differential equation 

3 

with B as defined in Eq. 1|12|) (with h = 0). This equa- 
tion, which can easily be derived from Eq. ©, was used 
to project from a solution found at /3 to a good initial con- 
figuration at (3 + A/3. Using this method we were able to 
track the solution of Eq. © over the temperature range 
(3 = 0.1 . . . 100. 

We have found it useful to split Eq. ||SJ at h = into 
two parts corresponding to the eigenvalues that are going 
to zero as (3 — > oo and the rest, 

, E M! + ? Mf. (17) 

A„ going to X n staying finite 

In a finite system at finite temperature, all eigenvalues 
are naturally non-zero, such that this equation is well- 
defined. However, for (3 > f3 c the eigenvalues which be- 
come exactly at T = decrease with system size as TV ~ 6 
at finite temperature, where b is an unknown positive 
exponent, and so are equal to in the thermodynamic 
limit throughout the whole low temperature phase. In 
this situation, the second sum in Eq. H17fl is equal to the 
diagonal elements of the Moore-Penrose inverse (see, e.g., 
[l2f ) of A which corresponds to the physical susceptibil- 
ity xa = d(sf)/dhf = P((sfsf) - (s?)(s«)). For f3 > f3 c 
Eq. 117|) can therefore be written as 

P = P(s?) 2 +Xu, (18) 

while for (3 < fl c , xu = Xu = P- I n Fig. 01 we have 
plotted x, defined as Xu averaged over sites, and its vari- 
ance Var(xii) as a function of (3. The plot shows that 
above the transition temperature x is equal to (3, whereas 
below the transition temperature \ remains essentially 
frozen at a value of (3 C . The site-to-site variance of x 
decreases with system size, roughly following a power 
law Var(x ri ) ~ iV -1 / 3 (data not shown here). Eq. (|T%|) 
then implies that the Edwards- Anderson order parame- 
ter q = (sf) 2 is equal to 1 - (3 c //3 = 1 - T/T c , i.e. the 
frozen components of the spin are those associated with 
the null eigenvalues of the matrix A. 

We believe that the approach to spin glasses via a 
1/m expansion method, while it is not simple to carry 
out, seems to be the only calculational technique which 
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FIG. 3: The site-averaged susceptibility \ (upper curves, left 
axis) and its variance (lower curves, right axis). The solid line 
corresponds to \ — P f° r comparison. 



can avoid the problems associated with the unwanted [2| 
replica symmetry breaking of the loop expansion. The 
point about which the expansion takes place, viz the 
large-m limit, is also interesting in its own right as the un- 
usual phase transition mechanism found in the SK limit 
would be expected to carry over to finite dimensions. The 
work by Viana \\?\ , who found that in the large-m limit 
the upper critical dimension and the lower critical di- 
mension are both equal to 8, indicates that below the 
upper critical dimension there is no Edwards-Anderson 
order but instead perhaps chiral order as suggested by 
Kawamura EH (but contested by 0). Our approach 
might in the future be developed into a tool to study this 
controversy from a new perspective. 
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